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Bibliografia 


— Simulação numérica da velocidade instantânea em escoamento 
turbulento, Ordenada — V; abcissa — tempo. 

À curva corresponde a uma solução do sistema de equações 
d u; du; d*u, 
“-+du, Sa + E E (x) (i=1,2,3) 


z 
dt X, 0x, 


obtido por J. D. Domingos, 

Representação gráfica da velocidade num pento feito através da 
impressora do sistema IBM 360/44 do Centro da lIniversidade Técnica 
de Lisboa usando uma subrotina desenvolvidas por 5. D. Antunes 
(trabalho Integrado num projecto de Investigação do 1. A.C. em curso 
no Núcleo de Estudos de Engenharia Mecânica). 
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POLUME XXXIII PUBLICAÇÃO MENSAL 


TURBULENCE MODELS AND THEIR EXPERIMENTAL 
VERIFICATION 


INTRODUCTION 


O. As an outcome of research work in turbulent 
flows started at Instituto Superior Técnico, 
Department of Mechanical Engineering, Division 
of Applied Thermodynamics (sponsored by Ins- 
tituto de Alta Cultura and Fomento da Acção 
Educativa do M. E. N.) with the close cooperation 
of Imperial College, (University of London, 
Department of Mechanical Engineering, Heat 
Transfer Section) and especially of Dr. J. H. 
Whitelaw, an extension course on «Heat and 
Mass Transfer in Turbulent Bondary-layers» 
was given at Instituto Superior Técnico in July 
1969, by Dr. J. H. Whitelaw (1.C.) and J. J. D. 
Domingos. 

This course, centered on a presentation of the 
Spalding-Patankar method of computation of 
turbulent boundary layers together with compu- 
tational results of Engineering importance, was 
a modified form of a post-experience course 
given at Imperial College in December 1968 and 
December 1969. The main contents of this course 
are now widely known through ref. (2)*. 

Due to the success of the course of July 1969 
at 1.5. T., a second course was evolved for 1970 
specifically to satisfy interests of the Lisbon 
audience, and to require a greater contribution 
from the local research group. Because much of 
the content of this course is not available in 
published form and because it was found 
desirable for future courses to have it available, 
the more inaccessible material is to be published 
in «Técnica» having in mind the wider audience 


* References collected at the end of part II. 
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which can benefit from the material and at the 
same time, became acquainted with work presently 
in progress at 1.5. T, 


0.1. The papers which follow (Parts 1, II, III, 
IV) closely resemble the lecture notes specially 
prepared for the course. The presentation could 
undoubtedly be reworked to improve clarity and 
comprehensiveness but this would imply delaying 
publication for am even longer time than has 
already elapsed. Neverthless, we expect the 
material to be useful to potential readers who 
may consider it as introduction to a increasingly 
important subject. If, as we hope, new courses 
follow the two already given, we confidently 
expect this publication to be a forerunner for 
more ambitious ones. 

0.2. As originally given, the course included 
lectures and experimental demonstrations which 
included the preparation of hot wire probes and 
manipulation of the equipment by participants. 

The lectures where given by J. J. D. Domin- 
gos (I.S.T.), J. H. Whitelaw (1.C.), F. Durst 
(1. C.) and S. D. Antunes (I.S.T.) and the 
experimental sections were conducted by F. 
Durst, S. D. Antunes and M. N. Mina. 


1. Once the «how» is explained, if the patience 
of the reader is not already exausted, the «why» 
of the contents is possibly worth further comment. 

The «why» comes from what we believe must 
be known by newcomers to the field, and what 
we think should be remembered by experienced 
readers. 

For the mathematician, turbulent flow in incom- 
pressible fluids only require the solution of 
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a set of four quasi-linear parabolic equations 
known for more than one hundred years but 
which no-one has yet been able to solve without 
drastic simplifications. And from the time of 
Poisson, Navier and Stokes the best mathema- 
ticians in each generation have abandoned the 
problem as untractable and apparently unsolvable. 

For the physicist, turbulent flow is the only 
great unsolved problem of classical physics... 

The engineer, however, is obliged to consider 
fluid flow in many aspects of daily routine. He 
designes aeroplanes, tall structures taking account 
of wind loads, turbomachinery, dams, space 
vehicles, artificial heart-lung machines, heat- 
-exchangers, combustion equipment, etc. with 
the constant challenge to produce better and 
cheaper designs. He can not abandon these 
problems. In the good old days, crude experi- 
ments and empiricism were enough. Later, 
dimensional analysis helped. Today is still the 
main tool... but no longer sufices. 

Yet, the potencial means of solving turbulent 
flow problems is still the solution of a set of 
four partial differential equations, or five, or six, 
or any number, if we add compressibility, heat 
transfer, chemical reaction ... for which a comp- 
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lete proof of physical validity is still to be 
produced. 

If one remember this, and also that such 
famous names as Landau, Onsager, Sommerfeld, 
G.1. Taylor, Kolmogorov and Kampé Forrié (some 
of them Nobel prize winners) abandoned the 
problem for more fruitfull areas of fundamental 
research, one must refrain from a frontal attack 
on some crude simplifications currently in use. 

There is only aestectic value, if any, in 
proposing solutions in terms of constants, or 
integrals etc, whose translation to nmumerical 
answers is beyond human reach; they do not 
give any fundamental insight into the physical 
mechanisms involved, besides the personal 
believe of the proponent in his own physical 
intuition. 

Summing up, we would like to escape in the 
20th century what Sir Cyril Hinshelwood said 
about 19th century: «Fluid dynamicists were 
divided into hydraulic enginners who observed 
what could not be explained, and mathematicians 
who explained things that could not be observed». 


J. ]. D. Domingos 
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|I— THE FUNDAMENTAL EQUATIONS OF TURBULENT 
INCOMPRESSIBLE FLOWS 


RESUMO 


Dedusem-se as equações fundamentais para o estudo 
dos escoamentos turbulentos incom pressiveis e discutem se 
as hipóteses fundamentais e suas limitações. 

Apresenta-se um resumo sucinto da posição actual da 
teoria matemática das equações de Navier-Stokes e do 
seu tratamento numérico. 

Discute-se em pormenor a equação de energia cinética 
da turbulência e a sua correcta dedução, tendo em conta 
o papel central que a mesma desempenha na formulação 
de modelos de turbulência. 

Todo o trabalho é orientado sob o ponto de vista da 
engenharia, no sentido do cálculo do escoamento médio. 


o. BASIC ASSUMPTIONS AND RELATIONS 


Except when explicitly stated, we assume : 

— homogeneous, incompressible and Newtonian 
fluids 

— isothermal flow 

— cartesian frame, fixed to the Laboratory. 

The substantive derivative (or derivative 


following the motion) is expressed by e 
t 


LM 
dx; 


Rd 
dt dt 
where, as everywhere, a repeated index means 
summation. 

The Reynolds transport theorem is expressed by: 


1) + V 


Ê d 6 d oNV, 
O Sa D+ To 
and mass conservation by: 
dV o V, a é e 
3 = —— sm ecause o 
) dx, dx, 
incompressibility) 
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Professor of Mechanical Engineering 
Division of Applied Thermodynamies 
Instituto Superior Técnico. 
LISBON 


SYNOPSIS 


The fundamental equations for the study of turbulent 
incompressible flows are derived and its underlying 
assumptions and limitations discussed. Present know- 
ledge regarding mathematical and numerical aspects of 
the Navier-Stokes equations are briefly presented. Special 
attention is given to a correct derivation and discussion 
of the turbulent kinetic energy equation having in mind 
its central position in the development of turbulence 
models for prediction purposes. 

All the paper is oriented towards engineering predic- 
tion procedures of mean flow behavionr. 


0.1. General form of a balance equation : 


dó) 
4) de dx, Jo); — Og 
Usually : 
Qg = G—s 


with G >> 0 — production of é (positive source) 
* > 0 — dissipation of d (negative source) 
if 
OQ; = G—:=0, ó is a conserved property 
For a volume element fixed in space (Euler 
description), we have, from 4) and 1) 


oq ) 
6) E rd dx, (Y, é + Jg) — Og 
where use has been made of the identity 
oú = d d V; 
7) ho 0x, se dx, O Nj) om dx; 
o 5) 
or 8) V Ox, =— 0x; (6 V;) 


because of 3). 


The term & V, is the convective (or advective) 
flux of é, and is purely due to the fluid motion, 
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while Jó; is a flux of é independent of that mass 


motion. 
From 4) and 1) 


dy ou o 
Mi Ox 


fo) : 
) ot 


0 
Ox, Ug)—e (G—s) 


1.0. The Navier-Stokes equations 


For a volume element 


dV, OZ, 
10) p =-——— +, 
dt dx, 


The terms which appear have the following 
physical interpretation: 


d V, 
- F “== mass times aceleration in the i direction 
t 


Z, — stress tensor: stresses applied for unit surface 
of the volume element (tensorial density) 


d ; . : 
uma cv" — resultant force on direction i of 


) 
Ux, 


the stresses on the surface 


: F -—i component of the applied forces for 
unit mass 


Equation 10) can also be expressed as a 
momentum balance. 


From 2) 


dV, d Ê ) 
Ro Ade AS, 


from 10) with 8) and putting 6=V, 
7, 


EO Cê md io 


which expresses the local (Euler) balance of 
momentum. 

Equation (10) or 11) are indiferently used as 
starting equations. 


7, 
dx; LeViV,+ TpJ+er, 


1.1. For an incompressible Newtonian fluid 
the stress tensor has the form: 


p+ q 4 d Vi RT 7, V, 

Tu P*O, Pa dg a e 

where P* is the hydrostatic pressure. 
From 11) and 12) 
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) o é 
de tevp= = O) Pesgterv VT 
13) ot ' 0x; 
( oV, 4 o V, LE 
Lãs q 0X; dx; ) ii 


where the terms have the following interpretation: 


É 

E” (2 V,) — rate of accumulation of momentum 

in the i direction in the volume 
element 


P+5, + V,V,= reversible flux of momentum 


OV,  OV 
a ( + 


) == momentum flux due to 
d Xj 0X; 


viscous forces 


For the incompressible fluid with constant 
viscosity 


d oNV, d dV, 
——— === = | — 
dx; dx, dx; dx; 

d oNV, 
0X; dx; 


oV, 
because = O 


x. 
So, we can write equation 13) as 


uj Eiji DAS pr 
dt dx; 1) 0X; 
or 
oV, 9V, 
15) -—— + j — 
t 0x; 


92V 
1 É. sa PR 


0 0x, d x 


Equation 14) or 15) with 


d Vv 
16) = 0 
x, 


are the basic equations for the motion of the 
Newtonian incompressible fluid, or, in short, the 
N.S. equations. 
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2. THE KINETIC ENERGY EQUATION 


Some of the terms which appear in the N,5. 
equations are better understood when the kinetic 
energy balance is considered 


For simplicity we assume F, = O 


The kinetic energy per unit volume is 


and its change in time is 


) IV 
dE lê [= PRN = eV, Mn: 5 
dt ot VA dt 
V, 
Taking from 15) and writing 
o V, (*) 
17) p= — 
dx, 
we get 
18) R 1 VV) 
ot [á i 
* d € 
amy são pos 
d X; 0X, dx; 


which is the balance equation for the kinetic 
energy. 


Taking into account the identities: 
V, r E ? la À 
dx; o x 


(*) The general relation is 


=| —— + =) , as is wellknown. 
! x. 


However on the basic equation Ya only enters through 


J 
divergence (= e T) which, for the incompressible fluid 
N j / 


d V 
with constant viscosity reduce" to ju 1 as alreadyshown. 
d x 
j 
DV, 
with the 


dx. 
j 


advantage of avoiding terms that would havet be simpli- 
fied afterwards. 


This allows one to proceed as if tj =b 
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) V, D) 1 


VV— = V,=-— (—— V,V) 
dE 0x, " dx à *% 
O V 
cf E EM | 
0X; 2 2 dx; 
and, for the incompressible fluid, ——— = 0 
x 
j 
dV d 2 
VV, — = (= 
dx, 0x, 2 
d P* 
A (PV) 
OX, d x 


The kinetic energy balance equation takes the 
form : 


1 
j 
p+* y oV, 
Eme Cs) =st 
ca g ) E ii dx, 


where the terms have the interpretation : 


ie (: e V, vi) = change of local kinetic 
energy per unit volume 
and time. 
1 pr Rd 
Vl NE E) = energy (kinetic plus pres- 
i 2 > sá 2 Mia 
* “ sure) flux density in direc- 
tion j due to the mass 
flux in the same direction. 
Is a convection (or advec- 
tion) of energy. 
V; Ty =energy flux density in 
direction j due to the 
viscous forces. 


The first term on the r.h.s. is the divergence 
of the flux of the above quantities, and as such 
represents the difference among outflowing and 
inflowing quantities. The second term on the 
r.h.s. is a source term (negative). So 


oV, 
Ti; —— — is the dissipation of Kinetic energy 


for unit volume and time. 


As we see, the viscous stresses give rise to 
two effects: one is an energy flux (V, 7,) the 
other a dissipation. 

The energy flux due to viscous stresses results 
from the fact that the slowers (faster) portions 
of the fluid in contact with the fasters (slowers) are 
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accelerated (desaccelerated) due to the existence of 
viscosity. The result is a transfer of momentum 
which involves, obviously, the simultaneous 
transfer of kinetic energy. 


dV 
The term 7, E Ti as been called a dissipation 
ai 


due to its non-divergence form in 19). However, 
to be a dissipation it must be intrinsecally positive. 
Its positiveness results clearly from: 

d V d V, d V, ) 

im Es pm - = u 


Ta, - 
H) dx, X ox 


Besides, if we integrate 19) on the whole fluid 
and if we assume it to be bounded and V. = O 
on the bounding surface, we get, using the 
divergence theorem: 


fear o (88) 
x) fdc dx, dv 
v v 


' À V. 2 
which shows that | Ta ) is clearly a dis- 
( 4 


sipation. By the second law of thermodynamics 
E >O. 

For future reference we now write 19) taking 
into account 


d 


) l 2 
ES ss O ê ES A Aee e 
dt 2 dx, 
V; o P* v o 3 fi) 
“ P 9x Ea ” em 


2.1. Kinetic energy balance in a laminar 
boundary-layer 


To clarify the behaviour of the several terms 
in the kinetic energy equation we will consider 
now as a typical flow a laminar boundary layer. 
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Using the typical symbols of the boundary- 
-layer literature we call 


X| =X V,=U 
X9=Y and VW=V 
Ka a Va == VA 


Boundary-layer flows are defined as two dimen- 

sional flows where a dominant direction exist 

such that variations of any flow property are 
most significant in the cross direction. For plane 

flows, such «definition» allows (by an order of 
magnitude argument) the N.5. equations to sim- 

plify in plane steady flows to: 


) d d P* > m 
dU +V U E P | e 
dx dy p d x dy 
d ) 
e 4 OU =0 
d X dy 


the flow beeing on the x direction and a wall 
exists for 


x > O y=0 
fU=Uco = constfory=<<,x>0 
aP 
dx 


and the equation for the kinetic energy reduces 
accordingly to 


1 9 U? 1 dUÍ 
— UU + —mV— sa 
2 dx 2 dy 
v o“ U? “OU ê 
MS 2 + v dA == 0 
2 dy dy 


Which we can make nondimensional multiplying 
by —— 
od mal 


For an initial profile U==U for x=0 we get 
the representation of fig. 1 
| 
| 
| 


The existence of the wall and of a non sliping 
condition implies U=0 at y=-0 from which 
results a velocity gradient and a viscous loss, 
across the flow. However, the non-dimensional 
viscous dissipation is 
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KINETIC ENERGY BALANCE IN A LAMINAR 
BOUNDARY LAYER 


12 pps E 


[| 

mo a 
am 
nm 


A 
E 
" 


=12 


| 
Y/px 
1) Ro 1 (U JU ourA 
Voo 3 2 d x dy 
+ d2 U? 
RR ( Dia 
Uco3 " dy? | 
e. a 
) U 3 v(5y 
FIG. 1 
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and the kinetic energy loss at each point is 


which do not balance everywhere the dis- 
sipation. Summing the difference across the layer : 


2 
f pm sc Er dy = 0, so, a global balance 
O U. 2 dy 
exists 
between the total loss of kinetic energy and the 
dissipation, the difusion transfer (curell 2) just 
taking into account the local unbalance, transfe- 
ring the Kinetic energy from a point to the other. 
As fig. 1 shows, the dissipation is maximum 
at the wall, where the kinetic energy is smallest. 
The dissipated energy came from the outer por- 
tions of the boundary layer where the loss of 
kinetic energy was in excess of the dissipation. 


3. SOME MATHEMATICAL RESULTS CON- 
CERNING THE NAVIER-STOKES EQUA- 
TION 


3.1. The fundamental equations 15) and 16) 
are a coupled set of four equations with four 
unknowns: V,, Vo, V;, P*. 


Mathematically, we have a sistem of quasi- 
linear partial differential equations of the para- 
bolic type. 

These fundamental equations are by themselves 
far from beeing accepted as truly general and 
fundamental though many people seems to 
paraphase Langrange (changing Euler to N.5S.) 
when he wrote in 1788: 

«One owes to Euler the first general formulas 
for fluid motion... presented in the simple and 
luminous notation of partial differences... By 
this discovery, all fluid mechanics was reduced 
to a single point of analysis, and if the equations 
involved were integrable, one could determine 
completely, in all cases, the motion of a fluid 
moved by any forces...» 

Indeed, before general solutions were found 
the Euler equations shown to be at most an 
approximation. Our hope have, after, been 
transfered to the N.S. equations. We have not 
yet found any general (which mean three dimen- 
sional, time dependent) solution, but we know 
already that they can not be general. They are 
better than the Euler equation, but we don't 


know yet to what extend. However, exact 
analytical results are of outstanding value even 
when they concern only incompressible fluids 
which are a mathematical abstraction without 
thermodynamic validity. 

For incompressible fluids existence and unicity 
of solution has been prooved for two-dimensional 
flows. It has also been proved that a steady 
state solution exists which is the assimptotic 
limit of a time dependent solution when time 
goes to =, 

For three-dimensional flows no such proof is 
known. Under certain restrictive assumptions, 
unicity of solution exist, and if the flow is 
completely bounded by a solid wall with non 
slip conditions, using methods of functional 
analysis it has been shown that the velocity tends 
to zero with time : This is a trivial conclusion from 
a physical point of view as the equation for the 
kinetic energy shows. 

The rigorous proof that N.S. equations admit 
turbulent solutions has not yet emerged, neither 
the conclusion that they do not admit, thougth 
some results point to the first possibility. At 
least, a very interesting conclusion exists already: 
two dimensional flows do not admit turbulent 
solutions which is in complete agreement with 
the physical experience. 


3.2. The almost untractable mathematical pro- 
blem, with the tools which have been used, does 
not prevent, however, some interesting conclu- 
sions if some additional assumptions are made. 
These are the object of the following chapters. 
Before, a comment must be done regarding 
numerical methods: even with present day com- 
puters the solution of the full set of time-depen- 
dent equations is impossible — the restrictions 
put on the stability and accuracy of presently 
known mumerical methods grow so sharply with 
the Reynolds number that the memory-speed 
requeriments of the computer go beyond the 
foreseen future. A crude, lower bound estimate 
shows these better than words; 

Let dU, dx,», be the caracteristic velocity 
variation between two nodes, 4x his distance, v 
the kinematic viscosity. 

Numerical stability requires the cell Reynolds 
number : 


AU Ax 


v 


<=1 
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If a flow is divided in 20 nodes in each direc- 
tion, the number of nodes will be 20>< 20 x 
>< 20 = 8.10 . But this number is at most accep- 
table for completely laminar flow on the Stokes 
or Oseen region. To go to an acceptable preci- 
sion an order of magnitude at least is needed, 
i, e. 200 ><200 >< 200 = 8.10 nodes. This is on 
the order of 10" greater than possible with pre- 
sent day computers. 

This again shows that one must turn to cruder 
methods if at least an engineering solution is 
though which is not the crudest. 

The way is to regularise the solutions with 
physically based guesses. 


3.3, Turbulent Flow 


As is known since the basic experiments of 
O. Reynolds, above a certain value of the Rey- 
nolds number the flow becames unstable. The 
exact upper bound of this critical Reynolds num- 
ber is unknown because it depends on the expe- 
rimental care, surface conditions, inicial profile 
(in a tube, for exemple, values between 2100 and 
70.000 have been found). The lower bound is 
much better known (i.e, for a tube for Re < 2000 
no turbulent flow exists). 

Again the question arises as to the validity 
of the N.S. equations. The stability theory of 
these equations shows the existence of an insta- 
bility of the flow in the sense that there are 
conditions in which a small disturbance is ampli- 
fied, not damped. The most success full approach 
to this stability problem has been the «linear 
stability theory» where the N. S. equations are 
linearized and on them the evolution in time of 
a small perturbation considered. However, the 
growth of the perturbation can not be followed 
in this theory by the basic reason that when it 
grows the problem ceasses to be a linear one. 

The non-linear theory, which would deal with 
the problem with not to drastic simplifications, 
is not yet known to have solved a simple real 
situation. 

Another approach which is slowly emerging 
is the so-called bifurcation theory. 

Bifurcation theory on its present form regar- 
ding the N. S. equations is much more con- 
ceptual than mathematical, but is worth mentio- 
ning because it clarifies, at least, some common 
concepts. 
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This theory assumes thatfor a certain R, < R, 
the first unstable mode of the linear theory will 
not grow exponentialy (as the linear theory would 
imply) with time but would tend asse ptpticaliy 
to a value proportional to (R — R« 3 

So, at R > Rc, two solutions of he N.S. are 
possible: one, the laminar solution, the other 
a periodic one of small amplitude, proportional 
to (R — Rç) das 

This situation will continue until a second cri- 
tical Reynolds number, Rcois reached, where from 


the laminar solution a second unstable mode 
grows. The process is repeated for each of the 
possible solutions. 

In this concept, there are not one but a gro- 
wing number of critical Reynolds numbers, with 
R.. In this way, for a Reynolds number well 
above the first critical one the number of allowable 
solutions of the N. S. equations is enourmous for 
the given set of boundary conditions, with the 
flow possibly switching from one to other con- 
tinuously. From the experimental (macroscopic) 
point of view, the flow appears as completely 
random ; this is the region of fully develloped 
turbulence. 

On the sequence of the exposition we will 
always assume fully developed turbulence and 
will describe the flow by its mean properties. 
Besides, we will assume the N. S. equations as 
still valid, and accordingly will introduce in them 
an averaging process. 

The averaging process must be considered 
more of an heuristic than a rigorous mathematical 
operation. 


we define: * 


* we Could have defined an «ensemble average» 
coined after the statistical mechanics concept. This would 
be pretensious and at the risk of inducing in the reader a 
sense of rigour which in fact does not exist on the physical 
problem. As a concept, it allows a correct definition of 
averaging. But that averaging can only be realized on a 
mental experience. So, on the scope of this course we 
discard it and use, instead, the experimentaly measurable 
mean quantity defined by the following relations. These are 
only rigorous on the limit of 4t = co Hcwever àt = o0 
would surpress all unsteady turbulent flows from the 
discussion. For A t finite, the approach is to consider it as 
a typical «time scale» of the flow under consideration. 
This, again, is a typical dubious definition of present day 
turbulence theory ... 
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dt 


u,=U,—-V, and x,=0 


thatz2 
cd] Prdt=P* 


ai 3 t/2 
p=P*— and p=0 
If we do the substitution 
V=U,+u, 
P*=P+p 


in the N. S. equations (and this is possible on 
mathematical grounds whaever the definition of 
U, or u,..., because we are only doubling the 
number of unknowns): 


o (U, U,) 9(U +-u) 
am CE o 4 eg) qo 
“as dx, 
J ó* 
aa (PA) 4 Êo] a) 

o o 0x, 0x; 
(U, + u)=0 

Ux 


integrating now between t and t + T and assu- 
ming the averaging concept as acceptable and 


t+at2 
taking account that -” and / are 
nes JS t=Ath 
comutative operations (assuming continuity) we 
get, 


9U oU, du; 

dt dx; dx; 
2 
d 
ed PE sy 

E OX dx; dx 

oU, 
= 0 

dx 


where the overbar, denotes time average, alway. 


Taking into account that 


du, 
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we can write 23) as: 
d U, d U, 1 dP 
— + 
b 


J X; 


24) + U, = — 


dt Jd X, 


d oU, So é 
dj em v——— — UU, 
0x, dx, 


which are similar to the equations for purely 
laminar flow except for the term uu. This 
new term is not known because we have 
introduced u, as an aditional unknown in the 
starting equations. This aditional unknown which 
transform the system on an undeterminated one 
is the «price» for the mathematical regularization 
got, and it translates the loss of information 
introduced by the averaging process. 

The term uu, =R, is an analogue of the 
viscous stress tensor and is called the Reynolds 
tensor (disregarding a constant factor). In fact, 
regarding the mean flow , R, is a momentum 
flux due to the random motions of the fluid, 
analogue of the randon motion of molecules 
which account for the laminar viscosity. This 
analogy is, however, an overview because R,; 
expresses a macroscopic motion, which must obey 
the N.S. equations. As a first aproximation, 
started by Boussinesg, we can go further with 
the analogy and postulate 


d U, dd U, 


E os 


J X, J X, 


thereby introducing the concept of a turbulent 
viscosity 1”) 


Even as a crude assumption, », has many 
advantages from an engineering point of view 
and provided its limitations are known, this 
approach is one of the most sucessfull presently, 
even when a crudest assumption is made 
regarding », as in Prandtl mixing-length theory, 
to be discussed later. 

However, before we simplify the problem for 
engineering sake we must further explore the 
averaging process on the basic equations in 
order to get a better insight as to postulate R,,. 

The next step is to introduce the averaging 
process in the kinetic energy equation. 


4 — THE TURBULENT KINETIC ENERGY 
EQUATION 


41 Introducing 
V=Uru 
Pt =P+p 


in the kinetic energy equation 20) and taking 


into sccuunt 


VM = (U, +u/=U +4u, 


VV V =(U +u)(U +u) 


Í 


= U, U,U, + uuu+2uuU, 


we get 


= Poda | E A A EP A 9 (UU) 9U, OU, 
— 0 À); -16U. |— — | — — ee —— +, — 
ot [2 o 4 dx 1 il 7) 2 Ê 0x; a d%, dx 
O VE ste ) E p 1 du, u; du, ou; 

e is Fai) + dx peão e" li dx, dx E + 


(*) uv, is always considered as a scalar quantity which is 
physically unacceptable. That would imply R, = Q if the 
mean flow as no gradient. More than that, R; =u Ê is 


proporcional on the lh.s. to the mean kinetic energy, 
d U, 


but ui = = 0 if the flow is incompressible. This 
d X; 
would imply that an incompressible fluid would not have 
turbulence. 
With the Boussinesq approximation and » a scalar, the 


assumption is, at best, a usefull approximation. 
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If we compare 24) with the starting equation 20) 
we note that the first group of terms in 24) 
corresponds exactly to the kinetic energy of a 
flow with velocity U,, and the second group has 
the same meaning, in the mean, for a flow of 
velocity u;. 
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However, when the two flows are added, due 
to the nonlinearity, a new group of terms arise: 


) 1 — 
dx, [aU + U, | 


This group is just a consequence of separating 
the instantaneous velocity field in two by 


V=U,+u, 
with the aditional condition u, = 0 


This means that 24) is exact as long as the 
N. S. equations are, and the assumptions 
regarding continuity of the fields are verified, 

If we wish to separate 24) in the kKinetic 
energy of the mean motion and of the turbulent 
field, conceptual consistency with the definition 
u,=0 implies : 

-— The kinetic energy of the mean flow must 
result from the basic equations of the mean flow. 

— The kinetic energy of the turbulence must 
result from the equations for u, (before averaging). 

— The sum of the kinetic energy of both flows 
must be given by 24). 


For the mean flow we have, from 23) 


9U, JU, oP 
25) ot xo. 0x; 
9“U, d 
oU, 
== 0 
0X; 


and by diference to.22) 


du, du, 9U, du, 
26) & EA fe a me (), —— Es 
dt x » 


j 
0X, d j d 
p) 
dp. ou, 0 
dx; d dx, 0x; = 0x, Eu; u; 
du, o 
dx, 


From 25) we have already seen that relative 
to a laminar flow we have the additional term 
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a (cu,u,) (the divergence of the Reynolds 


stress tensor). When we multiply 25) by U,, rela- 
tive to a laminar flow arises a new term in the 
equation of the kinetic energy due to the work 
done by the Reynolds stresses: 


27) U 


9 oU, 
= guu, U|—çcuu, — 
dx; 0x; 


If we now compare 26) with the equations for 
a laminar flow of velocity u, we see that the mo- 
mentum balance of the total low V, =U, + u, 
and the definition of mean flow gave rise to the 
new terms: 


The last one is not «time dependent» (on the 
scale of averaging) and is just the divergence 
of the Reynolds stress tensor with the sign 
changed. 

Multiplying 26) by u, and averaging the first 
two terms in 28) give rise to the energy terms: 


oU, du, 
29) cu, u, + eU u, 
d X, J X, 


The last one can be written 


à a 
2 2 


which is easily seen to be a convective transport 
1 
of 
2 
Coming back to 24: we see that the terms on 
the last line have already appeared in 27) and 29). 
——oU, 


However, the term — cu, u, = in the equa- 
dx 
j 


tion of the kinetic energy for the mean flow, and 


pu) by the mean flow. 


Í 


the term gu u, 


in the equation of kinetic 
dx 


energy of the u, flow do not appear in equation 
24) of the kinetic energy of the total flow. These 
terms cancel out when both the kinetic energies 
are added and this clearly explains its meaning : 
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— they represent an energy transfer betwen the 
mean flow and the u, flow. 

Because de u, flow represents the turbulent 
field, we call the equation of its kinetic energy 
the turbulent kinetic energy equation. 


T ae 
Putting E uu, = q 


we can write the equation as: 


J 2 dq 
ua gq + — 
dt dx, 
di “a PD | =b — 
dx, | ( Hd dx, | 
30) ms 0U “ du, y 
ci OU. U té 
| | d dx ' 
| J 
d q du, 
Except for the terms U, ra and—çu, Ux 


the equation is exactly the same as for 

the laminar flow, and as such the terms have a 

completely similar meaning. A more familiar 

interpretation, however, is possible if we remark 
da? 

that the aditional term U, ——— is just a conve- 


J 
xs 


ction by the mean flow of the scalar quantity q 


This allows writing the 1.h.s. of 30) as 


dq” 
dt 
On the r.h.s. the meaning of the new term 
—— dU 
pu, u, — which is a source term has already 


JxX 
been aiphisied tt is a source term due to the 
transfer of energy between the mean flow and 
the turbulent field. 
du 


2 
Í e . 

) is the viscous 
d X, 


dissipation which is responsable for the conversion 
of the kinetic energy to internal energy*. 


The last term: “( 


* There are some confusions regarding dissipation Hinze 


du 


2 
says (p. 66) that ( ) is not a dissipation which 


o X, 
can only be due to typographical error. 

qe. 
— ) is the only intrinsically positive 


/ du, 

In fact a — 
dx 

i 


term on the Kinetic energy equation, Besides, in the 
absence of mean flow, integrating on a closed volume 
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4.2. In the usual texts the turbulence kine- 
tic energy equation is derived without refe- 
rence to the total kinetic energy. However, we 
must stress the fact that the basic equations are 
the N. S. time dependent equations and the 
kinetic energy of the total flow. This will avoid 
any misinterpretation regarding the physical 
meaning of some terms like the dissipation. When 
the splitting between a mean and a intantaneous 
flow is done new terms arise. There are terms to 
which a physical meaning can be atributed. For 
instance, a term equivalent to the viscous stress 
tensor-which gives rise to the equivalent of a 
dissipation, which is not a true dissipation 
because the Reynolds stresses are not viscous 
stresses. The analogy, however, is a very helpe- 
full one provided we are aware of its limitations. 

Because the use of vectorial identities and the 
continuity equation allows the equations to be 
writen in diferent form, the exact form of each 
term will not always be the same neither the 
physical meaning which we give to them. 
These must be born in mind when comparing 
diferent texts of different authors. 

We will present now a common form in many 
texts of the turbulence kinetic energy equation, 
From equation 26) if we write 


1 E se 9 9 
k= > uu =" au, (q =ek') 


Í 2 ! 


dk psi eat dU, 
31) = E dt te MM 
dt a dx, 
dk 1 dp du, 
+ U = — —— 1. “Su. 
| dx, 2 * dx dx, 


where the interpretation of the terms is: 


d 3 o 
d X; 1 d X, 


32) U, 


which goes to o or is bounded by closed walls (where 
U, = 0) gives 


d a IR ; du, 2 
a d qu. a. d 
dt | Men i x | x ) Í 


W. 


du 
| 
IÊ ( ) where not a dissipation no other 
dx 
j 


terms could be and the kinetic energy of turbulence would 
be conserved for ever. 
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e 


o 
õ dx, 


7 kº 


2 dx, 


dk* 
dx, 


=U +U +U 


This terms represent the convection or advection 
of kº by the mean flow. 
All the terms can be experimentaly found 


Impossible to measure directly with presently 
known measuring techniques. 


du ; PET 
0) cus We ( so) 


ER [a | , 
33) uu, io ane SU, 0U, | + uu, neo: 
dx, u,u, 3x, + uu, Dx, 0 Xz 
oU, 
Fa DU, Bi 1 sa == 
U, U; x + u,u, 5x. º dx; 
oU, 74 mma “Ds 
fa de Frau 2x, | Us Us ja 


These terms aré source terms and represent the 
«production» of turbulence energy from the mean 
flow, as already explained. This means a trans- 
fer of energy from one flow to the other. 

If the mean flow has no gradients (U = const.) 
there is no transfer of energy to the turbulence. 
So, unless turbulence is diffused or convected 
from other regions of the flow, the kinetic energy 
of turbulence vanishes due to the dissipation. 

If the flow is two-dimensional 33) reduces to 
the four terms enclosed by the dashed square: 
these terms will be discusssed later in connection 
with elliptic flows. 

All the terms in 33) can be found experi- 
mentaly. 


34) — — uu, = - 


É Eae 
+ dx, (us k”) 
These terms have the meaning of a difusion 


of kºby the turbulent motion itself. They are 
very difficult to measure experimentaly 


35) — u 


Is also the equivalent of a difusion due to the 
coupling betwen pressure and velocity flutuations. 
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9 
: 3 (k?) is difusion of k due 
dx 


to its own gradient, like heat difuses due to a 
temperature gradient and termal conductivity or 
momentum due to viscosity and velocity gradient. 


The first term, » 


The second term: 


du, ? du, Nº du, 14 
Co Md 
0x, dx; OX, 


du,N? du, du, É d us é 
E) CJ CJ) 

dx dx; 0X, 0X; 
represents the viscous dissipation of turbulence, 
as already discussed. 


du, N? 
Only the terms a ) are measured with 
dx, 


not to much dificulty. 


q? 
2 
d x; 


Usually, the terms 32) and » (kKº) are 


grouped under the heading of «difusion terms» 
and its value acessed experimentaly by diference 
in order to close the kinetic energy balance. 
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NOMENCLATURE 
E — Total kinetic energy per unit volume 
F — External force per unit mass 


q” — Turbuleént kinetic energy per unit volume 
k — Turbulent kinetic energy per unit mass 


(du u,) 
a +" 


Jg — Flux density of 6 

Qu— Source of é per unit mass 

U, — Mean velocity component in direction i 

u; — Flutuating velocity component in direction i 
V;— Velocity component in direction i 

X, 1X; (i,j==1,2,3) orthogonal coordinates 
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*; — Total stress tensor 
é — Any physical property per unit mass 


u — Dynamic viscosity 
v — Kinematic viscosity 
o — Specific mass 


Obs: — An overbar (i.e u,ju,) means always 
a time average — Repeated indices means sum- 
mation on corresponding terms (Einstein con- 


NV, ] ON, 
vention; de V-— = a MN 
dx. j dx 
j 
k 3 9"U, 3 O ) 
Uj/= 2 UU, ES e, == 
i=1 0X; j=1 0X; 0X; 
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| — PREDICTION 


DU. D. C. 932.517.4 :632.526 


PROCEDURES FOR TWO — DIMENSIONAL, 


BOUNDARY LAYER FLOWS: 
EQUATIONS AND ASSUMPTIONS 


RESUMO 


Consideram-se e discutem-se as hipóteses e equações 
relevantes para métodos de cáleuto dos escoamentos 
médios bidimensionais do tipo camada limiteem regime 
turbulento. 


Formulam-se é discutem-se os tipos de equações que 
descrevem as curacteristicas de turbulência, desde O 
modelo simples que não utiliza equações diferenciais, 
até aos mais complexos que envolvem sistemas de equa- 
ções simultâneas. As vantagens relativas dos vários tra- 
tamentos são apresentadas. A justificação [física para as 


hipóteses inerentes em cada conjunto de equações é indi- 
cada, 


1. INTRODUCTION 


The subject of the paper is clearly stated 
in the title. Prediction procedures of varying 
degrees of complexity will be synthesised using 
forms of the equations derived in the earlier 
paper of J.J. D. Domingos (1), together with 
physical reasoning. It is necessary to stress the 
words «prediction procedures». The subject of 
turbulence models is interesting for its own sake 
because it taxes our mental ingenuity and mathe- 
matical skill; it is useful in so far as it enables 
us to predict important flow properties which 
would otherwise have to be measured. In almost 
all cases measurements are expensive in time and 
money; the purpose of a prediction procedure 
is to replace expensive measurements by cheap, 
but accurate, predictions. 

Attention must also be drawn to the restrictions 
mentioned in the title. The reason for the res- 
triction of hoo-dimensional flow was mentioned 
in ref. 1; the solution of the three-dimensional 
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SYNOPSIS 


The equations and assumptions relevant to prediction 
procedures for turbulent, two-dimenstonal, boundary- 
layer flows are considered and discussed, Sets of 
equations, from the simple case where no differential 
equations ave used to describe the turbulence characteris- 
tics, to the more complex cases which require several 
differential equations, are formulated avd the relative 
advantages of each set is made clear. The physical 
justification for the various assumptions in each set of 
equations is indicated. 


equations requires more storage than is available 
in the present generation of computers. Boundary 
layer flows are considered because the equations 
which govern such flows are a reduced form 
of those which govern recirculating flows. Thus, 
the eguations discussed in this paper will be the 
parabolic, boundary layer equations rather than 
the elliptic steady Navier-Stokes equations. 
It should be remembered, however, that a very 
large class of How situations is governed by 
boundary-layer equations; the classe includes 
jets, wakes, mixing regions, wall wakes and wall 
jets, pipe flow, diffuser flow, etc. The subject 
of elliptic flows will be considered by J.J. D. 
Domingos in a later paper. 

In the subsequent sections, soluble sets of equa- 
tions will be synthesised to represent turbulent, 
two-dimensional boundary layer flow. Assum- 
ptions are necessary to obtain these soluble sets 
of equations and will be discussed in some 
detail. The relative advantages and disadvanta- 
ges of the various sets will be made clear. Spe- 
cific predictions will not be presented in this 


299 


paper but reference to sources of such predic- 
tions will be quoted. 

The five sets of equations to be considered 
here are noted in Table 1 and are classified 
according to the number of differential equations 
used to characterise the turbulent nature of the 
flow. By not specifying in Table 1 a number of 
equations for the turbulent flow greater than 3, 


it is implied that there is a law of diminishing 
return; this point will be clarified later in the 
paper. 

Of course, where heat and mass transfer are 
important, the equations for total enthalpy and 
species concentration must also be solved. The 
content of the present paper is restricted to the 
aerodynamic equations but a similar treatment 


TABLE 1 
| 
Number of diff-. 
erential equa- | ee endent cs 
variables of Auxiliary 
tions to chara- | es References Comment 
dincise: rim: | differential equations 
en equations 
x ? =pP|S0 20] > Simple and ble: well 
dv | 0Y imple ao re e: : 
tested: mixing length is 
H=cp5 Us 3,4 surprisingly general: no 
account of turbulence. 
etc 
— — dU 
u —uv=u 16 Influence of turbulence 
ay intensity is considered : no 
; k —uv=c,lk? | 18,23,24 influence of scale: not 
uma universal, 
T —uv= Ak 17 | 
k,k! H = cul ok 2 Signigicantly greater 
1 universality : significant 
- k,k/Iº py =Culapk ? increase in complexibility : 
k,S un =e,(2k)/S| 19,21,22 | different zero locations of 
| | uv and dU/dy not predicted. 
k,w H=cpkw 4 24 
kkluv 26 uv and dU/)y can now become 
3 Es NONE ai zero at different locations: 
| increased complexity. 
MORE 1 4 equations 27 The greater the number of 
23 equations 28 equations, the greater the 
universality. The nature of 
the necessary assumptions does 
not warrant this degree of 
| complexity for boundary-layer 
| flows. 
TETE e = | E: US 
su TÉCNICA Nº 406 


is possible for the total enthalpy or species con- 
centration; the present absence of reliable mea- 
surements, e. g. for vt, tends to render such an 
approach unjustifiable. 

As a final comment in this introduction, it is 
important to realise that, for present design pur- 
poses, only eddy viscosity models (zero differen- 
tial equations for the turbulence) have been tes- 
ted with sufficient care and over a sufficiently 
wide range. For instance, the mixing-length 
hypothesis, together with the computational pro- 
cedure of D.B. Spalding and S. V. Patankar (2), 
has been tested over a very wide range of flow 
configurations (*) and its successes and failures are 
now well known: this cannot be said for the 
higher order turbulence models. It is clear, 
however, that models such as the mixing length 
have defects and higher order models are neces- 
sary to remove the defects. Thus, the design 
procedures of the future are likely to utilise two 
or threee equations for the turbulent flow. 


2. ZERO EQUATION TURBULENCE MODELS 


In this case, the problem is to solve the mean 
flow equation: 


(1) 
oU, 1 0P o“ U, ) 
EE A em e rs O 
dx, o dx dx, dx 
in its boundary-layer form, 
(2) 
) ) o — 
U CU yes 4 Ef U Jd tv 
dy dx o dx dy dy 
together with the continuity equation 
oU, 
FR =0. (3) 


For practical purposes, the velocity profile deve- 
lopment, U=U (x, y), and the value of the 
wall shear stress, are the dependent variables 
of greatest interest. For laminar flow the last 
term in equation 2, the Reynolds stress term, 
is zero and similar or numerical solutions can 
be found provided boundary conditions are spe- 
cified at 


(*) This procedure was the subject of an extension 
course entitled «Heat and Mass Transfer in Turbulent 
Boundary Layers» at Instituto Superior Técnicoin July 1969. 
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x=0 , y=0 and y—o, 


The equation is parabolic, as are almost all 
equations considered here, so marching solution 
procedures are possible. In practice the stream 
function or non-dimensional stream function is 
used to satisfy equation 3 and eliminate V from 
equation 2; the resulting equation remains 
parabolic. 

For high Reynolds number turbulent flow, 
the Reynolds stress term of equation 2 is very 
much larger than the viscous stress term. By 
analogy with the viscous stress term, it may be 
assumed to have the form: 


ee ) dU 
dy CO gy (à ) : (4) 


The set of equations comprising 2, 3 and 4 is 
soluble provided » is specified in terms of the 
dependent or independent variables. 

Many eddy-viscosity hypotheses have been 
proposed and used, in conjunction with numerical 
procedures and further assumptions for the nea”- 
-wall region, to predict mean-velocity profiles 
and wall-shear stress. Both G. Mellor (3) and 
A. M. O. Smith and T. Cebeci (4) who are 
concerned principally with boundary layers on 
adiabatic plane walls, employ an eddy-viscosity 
hypothesis which, in the outer region of the 
flow, has the form: 


é Ta Cs, Uq 


where 2, (x) is the displacement thickness of the 
boundary layer. 

This particular hypothesis works well for the 
range of flows for which it was intended ; it is 
inappropriate, however, to the large class of 
flows which are characterised by maxima in their 
mean velocity profile. Clearly, therefore, this 
hypothesis is restricted in its applicability. 

A more general eddy-viscosity hypothesis is 
that based on Prandtl's mixing length, 

i.e. 


“ | dU 
”, = |" | — 
dy | 
or 
— old ) 
— uv=| TR Rad (5) 
dy dy 


where the mixing length, |, must be specified. 
This hypothesis with | specified by 
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0.09 
0.435 YG 
(6) 


|= 0.435y , 0O<ys 


| = 0.09 y,; cam  < 
“Cêa + qaas da SF 

has been widely tested at Imperial College and 
was discussed in the Extension Course referred 
to on page 301. In that course, it was shown that 
successful predictions of important dependent 
variables were possible over a wide range of 
dependent variables and for many pratically 
important situations. Nevertheless, it was clear 
that the empirical specification of mixing length 
and the hypothesis itself were not universal: for 
example, the constant 009 did not permit the 
correct prediction of both radial and plane wall- 
-jet flows. 

The above examples of eddy-viscosity and 
mixing-length hypotheses were presented to make 
the point that both are valid over limited ranges 
of How —the mixing-length over a surprisingly 
wide and useful range — but that both are far 
from universal descriptions of the mean motion 
of a turbulent flow. Consider also the following 
more general statements of the inadequacy of 
such models : 


1. Shear stress depends upon local mean velocity. 
But do upstream and downstream properties 
influence u, u,? 

2. Shear stress is zero when d U/dy is zero. Is 
this in accord with experimental evidence? 

3. Flow properties are uninfluenced by the scale 
and intensity of the turbulence. Can this 
be true? 


It is not true that the local shear stress depends 
only upon the local mean velocity or its gradient' 
The upstream history effects in boundary layer 
flows are wellknown and have been commented 
upon in ref. 5. On the other hand, the nature of 
parabolic differential equations does not allow 
the downstream properties influence the boun- 
dary layer flow. 

In strongly asymmetric flow, it is known, that 
the zero shear stress and zero local mean velocity 
gradient do not occur at the same location. For 
instance, in wall-jet flows, it has been shown 
that the zero shear stress position is located 
closer to the wall than is the zero local mean 
velocity gradient (6, 7). Also, in plane ducts 
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where one wall is rough and the other smooth, 
the zero shear stress location lies closer to the 
smooth wall than does the zero mean-velocity- 
gradient location (8). Finally, and mostimportant, it 
is well-known that the scale and intensity of 
turbulence can significantly influence mean flow 
properties of turbulent boundary layers in certain 
situations. Examples which may be cited include 
the very substantial increase in heat transfer 
rate caused by a small increase in free-stream 
turbulence close to the stagnation point of a 
cylinder (9); the enhanced mixing of a mixing 
region caused by an increase in free-stream 
turbulence (10); and the lower transitional 
Reynolds numbers which are accompanied by the 
higher turbulence intensity levels (11). 

The arguments which led Prandtl to propose his 
mixing-lengtn hypothesis are well-documented 
(12) and need not be repeated here. They imply, 
however, that any dissipation of turbulent energy 
is instantaneously countered by a production of 
an equal amount. This, in turn, implies that 
diffusion and convection of turbulent energy are 
absent. But experiments have shown that this is 
only true in the logarithmic region of a wall 
boundary-layer flow (12). This represents a 
further defect of the mixing-length model. 

An interesting exercise (13) which may be 
performed in connection with eddy-viscosity 
hypothesis is to differentiate the general form 
of the eddy-viscosity hypothesis, 


o U, d U, | 
t- ! (7) 


— WU. UU. == v 
dx; dx; 


, YU l 


with respect to x, and multiply the result by U,. 

This will yield a shear stress equation which 
may be compared with that derivable from the 
Navier-Stokes equations without the eddy- 
-viscosity hypothesis. For the two equations to 
be equivalent, certain terms must be equal. The 
reader should consider for himself whether this 
circumstance is likely. 

The broad conclusion which may be drawn 
from the above discussion and from applications 
of eddy-viscosity hypotheses * is that although 
they have a limited range of applicability, the 


* The reader is referred to reference 15 for details of 
some further successes and failures of eddy-viscosity 
hypotheses, 
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range is sufficiently wide to justify their use in 
many cases. In particular, the mixing-length 
hvpothesis has been shown (14) to permit 
satisfactory predictions of practically important 
dependent variables over a wide range of flow 
situations. Eddy-viscosity hypotheses are easily 
understood, relatively simple to use and, since 
they require the solution of only one differential 
equation, comparatively economical of computer 
time. 


3. ONE-EQUATION TURBULENCE MODELS 


In devising a turbulence model which requires 
the solution of one differential equation for the tur- 
bulent flow, perhaps the simplest and most obvious 
approach is that of L.S.G. Kovasznay and V. N. 
Nee (16). In this case, the authors solved 
equations 2 and 3 together with equation 4 and 
a differential equation for v. The differential 
equation for »! was not derived from first princi- 
ples but was written down in the form given by 
equation 8: 


ps VE [ue jd 


0x dy OdyN Oy 
convection diffusion 
(8) 
dU 
+ (A cdr q (vi —v) 
production dissipation 


The form of this equation is similar to that of 
other equations which may be derived for the 
turbulent flow (as will be shown later) and this 
was the part justification for the inclusion of four 
terms. The coefficients A and B are hcpefully 
universal constants and L is a length scale which 
must be specified empirically, The equation 
implies that the eddy viscosity is convected, 
diffused, produced and dissipated in a manner, 
for example, similar to turbulent kinetic energy. 
The criticisms levied at eddy-viscosity models 
also apply to this one but, since there are more 
constants to adjust, it is likely that a slightly 
greater degree of universality may be attained. 

The more useful one-turbulence equation mo- 
dels are those which include the solution of the 
turbulent kinetic energy equation in addition to 
mean momentum and continuity. The two dimen- 
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sional form of the kinetic energy equation, 
without assumptions, may be written as: 

dU d 3 
TE ed + —— u, u, e 
dx; 0X, 


u, OP du, 
"a e e a SA 
e dx; dx, 


The inclusion of equation (9) permits the mean 
flow to be influenced by the turbulent energy 
but the deficiencies 1 and 2 of the eddy-viscosity 
model still apply. The terms in equation (9) have 
been discussed in some detail in ref. 1 and more 
will be said, regarding their measurement, in 
a subsequent paper. The two-dimensional, boun- 
dary-layer form of equation (9) may be written as: 


U mm + 2u,u, 


dx, 


(9) 


d d ————— OU — OU 
U—k+V—k+(uy—v'))— +tuv— 

dx dy ade dx dy a 

convection production 


o 2 " 
+ ykrpo- kt) fe=0 
dy dyº 


— 02-01. — 


diffusion dissipation 


(10) 


I should be clear that the set of equations 
comprising continuity, mean momentum and 
equation (10) cannot be solved without assump- 
tions. First, a link must be provided between 
uv and k. One such link is that used by P. 
Bradshaw (17) i.e. uv « k; another is the Prandtl- 


-Kolmogorov formulation used extensively at 
| 


! 


Imperial College, i.e. uv k? dU/dy, or » = 
l 


= €,! sok”. In addition, the diffusion and 
dissipation terms are intractible in their present 
form. For high turbulent Reynolds numbers, the 
viscous dissipation may reasonably be neglected 
but the remaining component of diffusion remains 
and includes both v and p about which our set 
of equations provides no information. Indeed, 
the fluctuating pressure represents a substantially 
unknown quantity since it has not been measured 
except at a wall. The term is, therefore, approxi- 


mated by: 
h 7 v, dk 
dy kt dy 
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There is little direct evidence of the validity 
of this assumption which implies the existence 
of a Prandtl number for turbulent energy and 
requires that the kinetic energy diffuses down 
its own gradient. On the other hand, other 
scalar quantities such as enthalpy or mass con- 
centration behave in this way so, in the absence 
of alternative information, it is a plausible assum- 
ption. Of course, the plausibility decreases as 
eddy size increases. With this assumption, we 
must continually ask the question «is 7, , a cons- 
tant or does it vary from one flow situation 
to another, or from one part of the flow to ano- 
ther?» Tests carried out at Imperial College sug- 
gest that Cut is constant over a very Wide range 
of flow situations. Further testing is, however, 
necessary to confirm the actual range over which 
it may be assumed constant. 

The dissipation term, €, must also be appro- 
ximated. It is defined as: 


ra 
E = ( ) : (11) 
OX; 


du, N? 
The nine components of E cannot be 


Xj 


measured and cannot be included, per se, in 
a one-equation model. It is assumed, therefore, 
that energy is transferred from larger to small 
eddies (a cascade effect) until there are no small 
eddies to which the energy may be transferred. 
The very small eddies are assumed to lose energy 
as a result of viscous stresses. So dissipation 
takes place at high wave numbers (i. e high 
values of 27f/U) and in a region of the energy 
spectrum where it can be argued that the tur- 
bulence is isotropic. With this 'assumption it 
can be shown that: 


iss C, k é 
(: (12) 
lp 


where C,, is a constant and 1, is a length scale 
characteristic of the size of the larger eddy. Thus 
the diffusion and dissipation terms equation (9) 
may be rewritten as: 


) o. U — () 
U pap lado 4 uy o 
dx dy dx dy 


+ 


eo 
— = 0 (13) 


v dk 
4 d ( t bd A 

dy Ns dy lp 

For boundary layer flows, the production by 
normal stresses can usually be neglected and, if 
this is so, the set of equations comprising 
continuity, mean momentum and turbulent kinetic 
energy is soluble provided that a turbulent 
viscosity law is specified. D.B. Spalding (18) 
strongly favours the Prandtl-Kolmogorov visco- 
sity law in the form: 


k* (14) 


This set of equations requires a specification 
of Tp te Coe lp: € and | and these must be found 
from experimentalinformation. Because dissipation 
and production are known to be equal in the 
logarithmic region close to the wall, the law of 
the wall provides a link between C, and C, ,i. e. 


Es 


Es U n (— de as 


(15) 


+ '/a u 
Ca (7w/2) É 


Thus, a choice of C,, implies a choice of C,, and 
only four parameters need to be specified. 

Bradshaw's approach favours a viscositylaw lin- 
king uv with k only; this allows the kineticenergy 
equation to be transformed to an equation with 
shear stress as dependent variable. In this shear 
stress equation, the diffusion term is treated 
in a different manner from that indicated above 
and the resulting set of equations turn out to be 
hyperbolic; they are then solved by the method 
of characterisitics. 

There is little doubt that the one-equation mo- 
dels are superior to zero equation models from 
the point of view of generality; the larger num- 
ber of assignable constants ensures this. Brad- 
shaw's method has been shown to allow very 
satisfactory predictions of mean velocity and 
wall shear stress over a very wide range of boun- 
dary layer flows; the proportionality between 
shear stress and kinetic energy ensures, howe- 
ver, that the procedure will not work, in its 
present form, for flows with velocity maxima. 
On the other hand, the procedure was not deve- 
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loped with this in mind. It is also worthy of note 
that the mnumerical solution procedure discussed 
in the Extension Course referred to on page 17, 
is suitable for solving parabolic differential equa- 
tions such as equation (13). Thus, this numeri- 
cal procedure, which has already proved very 
satisfactory for zero equation turbulence models, 
is also applicable to one and to higher order 
turbulence models, provided the equations remain 
parabolic. 


4. TWO-EQUATION TURBULENCE MODELS 


In one-equation turbulence models, it is neces- 
sary to specify 1, (and 1,) empirically. As was 
the case for the mixing length, this and univer- 
sality are incompatible. The logical step forward 
is to provide a differential equation with 1, as 
dependent variable and determine if this raises 
any new problems. 

Such an equation could be obtained in the 
same manner as L.S.G. Kovasznay and V. N. 
Nee obtained their equation for »,; that is, to 
write down an equation of the form: 

d | 


o dis 
a Vo iria 5y) + 8 D (16) 


where the form of the generation and dissipation 
terms, G and D, are guessed. A more satisfactory 
procedure is to attempt to obtain an equation 
from first principles. It may be remembered that, in 
the previous section, it was assumed that dissipa- 
tion of turbulent energy could be replaced by: 


Co (kk) 
lp 

This implies that if an equation for dissipation 
can be derived, it can be related to |, and k pro- 
vided local isotropy of the dissipating eddies is 
assumed. 

An equation for dissipation was derived by 
F. H. Harlow and P. I. Nakayama (19) based 
on earlier derivations by P. Y. Chou (20). The 
derivation can be set down as follows: 

The equation: 


| ——— u. — — 
dt ) 0x; | dx ) 0x, 
“du, 1 dp du; 
e =. e — + : 17 
ú; 0x, o dx d xi ( 
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may be obtained by subtracting the time mean 
components from the general Navier-Stokes 
equations. Together with continuity it implies: 


du, ia | du; oU, du, 
Fr dx, o) dx, iá dx 
1 dp du, RR 
Cop ox To fox Mu) (8) 


If equation (18) is differentiated with respect 


du 
to x, and multiplied by a and thenaver aged 
! 


the result is: 


J* U, du, 9U, ou, du, 
dx,dx, ! dx, dx, dx, dx; 
(a) (b) 
oU, du. du o) g É du. N2- 
a: ia: Mia 095 2 (5) | + 
dx, 9X, 0X, | dx, E ([m) | 
(c) (d) 
ou, du, du, J 1 du, | 
dX, UX, 0x; | dx É 0X; E 
(e) (£) 
1 du, Jp d u, u du, 
pi o! dx, dx, dai. Ux; UxS Ux, di 
(g) (h) 


1 B “du. dp 
= — — (— )+ 


dx, q dx, dx, 


| da ET 
j” | 1 (du, | | du; 
| q 
dé Td ER pe E -— 
Í 
1x 2 N OX, / , d X; 0X, 


The two-dimensional, boundary-layer form of 
equation (19) is 


(u “> + o E (128) de. 


od XxX d y d y 


(d) (a) 
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1 du, p 


e dx, Ox, 0x, 


du, du; du, ds 
— — — + wy 
0X, dx, dx; dy 


(e) (£) (g) 


(20) 


o | ma 


/ É 2 
with S = — f a) 


This equation has dissipation as dependent 
variable but is intractable in its present form. 
(It should be noted that, for simplicity, density 
has been assumed as unity). To make it tractable 
the following assumptions can be made (see refs. 
19, 20 and 21): 


— term (a) may be neglected from order of magni- 
tude considerations 


— terms (b) and (c), from a loose kinetic energy 
analogy, may be postulated in the form 


— terms (e), (f) and (g) are diffusion terms and, as 
usual, may be approximated by a gradient term, 


Be 45 


Ca dy 


Pa 


—the last term in equation (20) may be appro- 
ximated by 
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With these assumptions the equation may be 
written, for high Reynolds numbers, as: 


Equation (21) has the same form as equation 
(17) and may be solved together with continuity, 
mean momentum and turbulent kinetic energy 
if the eddy viscosity is postulated in the form: 


(22) 


This model requires the specification of 7,, 75, 
C,,C,and C,,i. e. one morethan in the one- 
-equation model of the previous section ; the need 
to specify the length scale has, however, been 
eliminated. This last point is important because 
whereas Prandtl numbers may reasonably be 
expected to be constant, length scales cannot. 

It is the writer's opinion that this model repre- 
sents the points of diminishing return of boun- 
dary layer flows with the exception of those flows 
where strong asymmetry is present; and even 
in such asymmetric cases, it may well not be 
worth solving additional equations. The model 
has been tested at Imperial College by W. P. 
Jones and B. E. Launder (22) and appears to be 
satisfactory over a wide range of flows. 

K must be recognised that having derived, with 
some assumptions, an equation for S or 1,,, alter- 
native equations can be derived. For example, 


the equation for S may be regarded as an equa- 
E 


tion for - - Equations can also be derived 


D 
for EL Where m and n can take on any 
value. D. B. Spalding and K. H. Ng (23) have 
preferred an equation for (k,,) and have used 
it to predict the data assembled for the Stanford 
Conference and also for supersonic boundary 
layers; in both cases the results are excellent. 

Also within the framework of the two-dimen- 
sional models, D. B. Spalding (24) has employed 
an equation for fluctuating vorticity in place of 
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the dissipation equation; that is, the dependent 
variable is w=(du/dy— od v/0x). They emp- 
loy the equation in the form: 


Ow Ow d [E =) + 


w 
2 
+ /0U 
+ €; dc aid (23) 
d 
the k-equation in the form, 
U dk V dk d po Ok/2: 
2 dy 2 dy dy NG dy ) 
(24) 


1 
te a pod + Cs kw* = 0, 
dy 


and the Prandt]/Kclmogorov hypothesis in the 
form : 


O 


He = Cuekw. (25) 

Of the models discussed in this section, the 
k —S and k —k, models have been subjected to 
the most thorough testing. Both have proved 
highly successful, although it remains to apply 
them to some important flow situations such 
as those with strong density gradient. These 
models do, however, seem very promising and 
it is likely that design procedures of the future 
will utilise two-equation models. 


5. THREE-EQUATION TURBULENCE MODELS 


Two reasons exist for considering the addition 
of a third turbulence equation. It will inevitably 
allow greater generality since it will involve more 
constants (or functions) which can be assigned 
values based on experimental information. The 
second reason is that, in strongly asymmetric 
turbulent flows, it is known that the surfaces of 
maximum velocity and zero shear stress will not 
coincide. Sometimes the distance between these 
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surfaces is negligibly small compared with the 
width of the boundary layer but on other 
occasions the distance is not small and must be 
considered if accurate predictions are to be 
achieved. An example of a flow situation where 
the non-coexistance of zero shear and maximum 
velocity is important is flow in an annulus where 
one wall is smooth and the other rough. In such 
cases, it is necessary to add an equation for the 
shear stress uv. An appropriate set of equations is : 


oU OdV 
dx E in 
9U dU dP d — 
E e ci Da 
E é V dk — “JU 
E dx "3 dy * gy ” 
7) um, Ok/2 5) 
LE »y ) + = 0 (26) 
ds oS s dU 
na e a 6 
+ (— ) e E a 
dy o dy/ 2 k o 
d uv d uv ») E. duv 
» dx de E as (= a 
0U Suv 


This set of equations requires a specification 
Of 0,4, 964, 9., C,, C, €, and C,. It does not 
require the specification of the length scale and 
does not use the Prandtl-Kolmogorov viscosity 
law. It is useful to note that when convection 
and diffusion of shear stress are mad equal, 
equation (26) implies that: 


i.e. the Prandtl-Kolmogorov eddy viscosity law. 
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The above three-equation model has been 
tested with some succes by K. Hanjalic and 
B. E. Launder (25) using the computational pro- 
cedure which was described in the previous course. 
Wall-iet data, two-dimensional plane rough-chan- 
nel data and rough-annulus data have all been 
well predicted. 


6. MAOY-EQUATION TURBULENCE MODELS 


It should now be clear that a system of equa- 
tions can be built up which contains a very large 
number of equations. The question must also be 
asked, however, is it worth doing? For boun- 
dary-layer flows, there is likely to be little merit 
in a system of more than three equation for the 
turbulent flow; and, for the vast majority of 
cases, two equations should suffice. Larger num- 
bers of equations have been proposed and it use- 
ful to point out how these are utilised and to 
mention some points in their favour. It is likely 
that a greater number than three turbulence 
equations will be required in recirculating flows, 
where the more complex elliptic, time-average 
Navier- Stokes equations must be solved for the 
mean flow and the description of the included 
in the various equations. Thus, in this brief 
section, the comments are of some relevance 
to the separated flow situations. 

In equation (13), the production of turbulent 
energy by the mean fllow working against nor- 
mal stresses must be neglected in a two-or three- 
equation model. This is because there is no spe- 
cific information regarding the three fluctuating 
velocity components. This is justifed on the 
grounds that, in the vast majority of boundary- 
layer flow situations, the normal-stress produc- 
tion is smal compared to the shear-stress pro- 
duction. However, in a region close to separa- 
tion, this is not true. In such cases, it may be 
desirable to replace the turbulence kinetic energy 
equation by three equtions, one each for the 
three components of fluctuating velocity ; length 
scales can be formulated for each of these three 
equations; shear stress equations with further 
length scale equations can also be derived. Of 
course, the number of constants increases at a rate 
of approximately twice the number of equations 
if the high Reynolds number assumption is made 
and may be even larger if this assumption is not 
made. Thus, the optimisation of the constants 
will present formidable problems. 
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7. CONCLUSIONS 


I is useful to review here the conclusions 
which have been drawn in the text. They may 
be stated as follows: 


1. Zero-equation turbulence models are currently 
in use for design purposes. They offer greater 
generality than correlations or integral proce- 
dures and have been subjected to considerable 
testing and evaluation. 


2. Two-equation turbulence models apper to 
offer the best compromise of generality and 
simplicity but have not yet been tested over 
a wide enough range of flow situations to 
allow them to be recommended for design 
purposes. It is likely that, in the future, such 
design procedures will be used. 


3. In strongly asymmetric turbulent flows, three- 
-equation turbulence models may be necessary. 


4. It will be necessary in all or any of the above 
turbulence models to embody them in computer 
programs which are economical of both storage 
and computing time. 


NOMENCLATURE 


C Constants 
| Length sclale 


e RR 

k Turbulent kinetic energy |- u+v Em | 
2 

Fluctuating pressure . 


P 
P Mean pressure 


Ju, * 
S  Dissipation of turbulent energy [= &) |, 


IT Temperature 

u Fluctuating velocity in x-direction 

U Time mean velocity in x-direction 

v Fluctuating velocity in y-direction 

V Time mean velocity in y-direction 

w Time mean velocity in z-direction [also = 
ass QEU Oda 


x 
| Orthogonal coordinates: x in stream-wise 
direction 


2, Displacement thickness 
o Density 


O 
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é Dynamic viscosity 

v Kinematic viscosity 

q Prandtl or Schmidt number 

t Fluctuating temperature 

f Frequency 
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4 HERÓIS DA LUBRIFICAÇÃO > 


Ben-Hur 


Toda a gente sabe quem foi Ben-Hur, o quadriguista 
mais famoso do seu tempo. 

O que poucos sabem é que Ben-Hur era um dos mais 
antigos chefes de manutenção da Companhia Grega de Car- 
ros de Corrida. E foi por ele ser homem dotado de muita 
experiência que conseguiu «limpar» o Grand Prix Anual de 
Quadrigas, no Estádio de Antióquia. 

Messala, o grande favorito, representava as cores 
romanas. Ás apostas eram de 8 para 1. 

Ben não passava de um desconhecido e ninguém se 
atrevia a investir o seu rico dinheirinho numa incógnita 
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Ben-Hur 


daquelas. Mas Ben-Hur teve uma súbita inspiração na ma- 
nhã do próprio dia da corrida. 

Embora sabendo que os seus cavalos não corriam 
tanto como os de Messala, que eram árabes e tinham a 
culassa rebaixada, ocorreu-lhe que uma lubrificação eficaz 
das rodas das quadrigas talvez fosse uma vantagem. 

O problema pedia resolução urgente, mas... onde é 
que lhe seria possível obter uma assistência técnica rápida? 

Consultou a lista de fornecedores de lubrificantes, deci- 
diu-se Ben, depois de um meticuloso exame, a pedir uma 
mãozinha à MOBIL. E porquê? Porque a MOBIL tinha um 
corpo muito completo de pessoal habilitado. Tudo levava, 
pois, a crer que lhe enviariam um Delegado a tempo. 

Ben-Hur tinha razão de assim pensar. 

O Delegado da MOBIL veio a correr, o caminho todo, 
com um frasco de «Mobilinicum XVII — Massa Lubrificante 
Especial para Eixos de Carros de Corrida». 

A corrida começou e o resultado foi o previsto (o pre- 
visto por Ben-Hur, que a multidão, essa, apostara em 
Messala...). Foi a lubrificação dos eixos do carro de Ben 
que venceu, incontestávelmente, a prova. Messala perdeu | 
por seis carros e, desgostoso, acabou os seus dias a conduzir 
carrinhos de mão... 

Moral para Messala: em Roma sê romano, em corri- 
das de carros sê MOBIL! 

Moral para os responsáveis pela manutenção nos nos- 
sos dias: não é só o custo que identifica um lubrificante. 
A assistência que funciona por detrás dele é decisiva — e 
de que maneira! 

A MOBIL, excepcionalmente bem provida de Delegados 
competentes e tecnicamente muito bem orientados, pode 
garantir aos clientes a assistência de que eles necessitam. 

Messala tinha toda a assistência do seu lado, mas 
essa não era «assistência» MOBIL... 

Os Delegados da MOBIL estão ao vosso dispor — 
com os melhores lubrificantes que podeis comprar e com a 
competência técnica de que precisais. 

Eles estão sempre onde são precisos! 


E exe Ea Ea VA DS id vas o ETA ves vd NÃ AC ESA Ae VA 
ENE LPA LYRA VAL PRANE LA VA PA LA VR IÇA EVA LAS 


NO, 3 da Sei Ie 


= Oil 


Loto PaSAMA 
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CONCEPÇÃO FÍSICA E INTERPRETAÇÃO ANALÍTICA 
DAS PERDAS DIELÉCTRICAS 


RESUMO 


Estudam-se as perdas por condução, devidas à polari- 
zação e causadas por descargas parciais, Para isso anali- 
sam-se as variações da condutividade deelectrica e da permi- 
tividade e calcula-se o esquema equivalente de um dieléctrico 
com cavidades, 


1. INTRODUÇÃO 

Em qualquer sistema eléctrico, por mais simples 
que seja, surgem problemas de isolamento. 
O conhecimento dos materiais isolantes aplicáveis 
a fim de se conseguir o conveniente funciona- 
mento é uma preocupação do engenheiro que 
projecta e que executa. 

No que se refere aos materiais dieléctricos 
é essencial a compreensão das perdas. Existem 
hoje teorias que dão em boa parte resposta a essas 
inquietantes interrogações, permitindo-se uma 
representação física racional, enquadrada nas 
modernas concepções da constituição da matéria, 
e uma tradução analítica fundamental na descrição 
dos conceitos físicos respectivos. 

Um trabalho de sistematização é em última 
análise o objectivo do presente artigo. Para isso 
estabelecem-se as expressões clássicas para o 
dieléctrico ideal, introduzindo-se a noção de permi- 
tividade (!) complexa para atender às perdas de 
polarização [1], e descrevem-se as ideias basilares 
da teoria de Debye que deu elevado impulso ao 
estudo dos dieléctricos [2]. Finalmente elabora-se 
um modelo aproximado para se deduzirem as 
perdas de ionização [3]. 


(') Usa-se o termo permitividade para denominar a 
grandeza : normalmente designada constante die- 
léctrica por não ser de facto um valor «constante» 
(depende da frequência e da temperatura em geral), 


TECNICA Nº 46 


por 


HERMÍNIO DUARTE-RAMOS 


Colaborador Científico do Schering - Institut (T.U, Hannover) 
Engenheiro Electrotécnico (1, S. T.) 
Bolseiro da U, L, (Angola) 

SYNOPSIS 


The losses of conduction, of polarizsation and by partial 
discharges are studied. For that we analyse the variations 
of the conductivity and the permitivity, and an equivalent 
diagram of the partial discharges in voids is caleulated, 


2. CONDENSADOR ELÉCTRICO DE VÁCUO 


Seja um condensador cujo dieléctrico é o vácuo, 
portanto constituído sem matéria. Quando a 
tensão aplicada aos eléctrodos do condensador 
for v a carga eléctrica acumulada na superfície 
dos eléctrodos é Q=C, v, onde €, é a capaci- 
dade electroestática do condensador, sendo a 
intensidade de corrente no circuito exterior 
i=dQ/dt ou i=cC,.dv/dt. 

Se a tensão for alternada sinusoidal v = V sin mt, 
onde V é o seu valor de pico, então i = | cos mt 
com o C V ou, em simbologia complexa, 
[= jo em V, apresentando-se uma desfasagem da 
corrente de 7/2 rad em avanço relativamente 
à tensão e portanto a potência activa, definida 
pelo valor médio do produto v.i, é nula dado 
que as grandezas v e i são ambas sinusoidais 
puras. 

Admitindo um condensador plano a capacidade 
geométrica C,, assim chamada por se caracterizar 
apenas pelas dimensões do condensador, será 

o = &» A/d onde £, é uma constante universal 
:,=8,85 pFm—!, A éa área da face dos eléc- 
trodos e d a distância entre eles. 


Define-se densidade de fluxo eléctrico ou deslocamento 
eléctrico pela carga eléctrica por unidade de super- 
fície D=q, ondeqg=Q/A sefor Q a carga 
total de cada placa do condensador (em rigor, 


pela lei de Gauss, O = bs D. dS onde S é qual- 


309 


quer superfície fechada que envolve a carga Q), 


— 


sendo vectorialmente D, =, E,, isto é, o deslo- 
camento eléctrico está em fase com o campo 
eléctrico. 

Sob o ponto de vista prático o ar comporta-se 
com aproximação suficiente, se não ocorrerem 
fenómenos de disrupção completa ou incompleta, 
como se se tratasse de vácuo: para o ar à 
pressão atmosférica a sua permitividade é 0,69/0% 
maior que o valor £, do vazio. 


3. DIELÉCTRICO PERFEITO 


Introduza-se no condensador um dielétrico de 
permitividade relativa : e de condutividade 
eléctrica nula (dieléctrico perfeito). Fisicamente a 
introdução deste material em substituição do vácuo 
provoca uma redução do campo eléctrico de :, ou 


seja E=E| /:, sendo D=-D,. Este facto resulta 
da polarização do dieléctrico sob o efeito do campo 
aplicado exterior que, orientando os dipolos 
elementares da sua estrutura na direcção do 


campo, provoca um campo de polarização E p,em 
—+> BEN 
oposição ao campo aplicado, tal que E=E, — 
— Ep = Et, 
Mas para a mesma tensão aplicada mantém-se 


o campo eléctrico (E=E,) pelo que deverá 
aumentar de c: a densidade do fluxo eléctrico, 


D=: D, donde D=:<, : E. O deslocamento 
eléctrico corresponde igualmente a uma densidade 
de fluxo eléctrico no vazio sobreposta a uma 
densidade de fluxó de polarização, chamada 


vector de polarização P=D—D,, donde D= 


= ES Pou D=< «E com P=+ (:—1).E, 
onde y=:—1 se chama susceptibilidade die- 


léctrica de modo que Poem Ec E. 

Daqui resulta que para se manter igual a 
intensidade do campo eléctrico com o vazio ou um 
dado material de permitividade : exige-se uma 
carga eléctrica : vezes maior no caso de existir 
matéria. Os dipolos de facto constituem longas 
cadeias de partículas, alinhadas devido à polari- 
zação, terminando com uma carga eléctrica posi- 
tiva junto do eléctrodo negativo e uma carga 
eléctrica negativa na extremidade próxima do eléc- 
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trodo positivo. Estas cargas eléctricas superficiais 
do dieléctrico, opostas às correspondentes cargas 
eléctricas dos eléctrodos, neutralizam parte da 
carga eléctrica das placas. Se for Q a carga eléc- 
trica total apenas uma parte Q/: contribui para 
neutralizar a tensão aplicada sendo a restante 
carga eléctrica Q (1 — 1/:) neutralizada pela 
polarização. 

Pela definição de deslocamento eléctrico é 


D = q e por conseguinte q dA = D. dS, onde 5 
é qualquer superfície que envolve A (lei de Gauss). 
A carga eléctrica relativa ao campo aplicado q/: 


no caso do vácuo conduz a (q/:). dA =, E. dS. 
A restante carga eléctrica define analogamente o 
vector de polarização P tal que q (1 — 1/:) dA = 
=— P. dS pelo que q dA =: E. dS + P. dS 
ou sejaD=+: E+P. 

O fenómeno da polarização equivale a um 
aumento da capacidade do valor :, isto é, C=- 
= C,, continuando a verificar-se as mesmas rela- 
ções para a carga eléctrica e corrente em função 
da tensão tendo em atenção o novo valor da per- 
mitividade, O == Cv, i=C. dv/dt, [= jo C V. 

Devido contudo ao fenómeno de polarização, 
mesmo num dieléctrico perfeito (de condutividade 
eléctrica nula), existem as chamadas perdas dieléc- 
tricas de polarização que se traduzem numa compo- 
nente óhmica da corrente. Este facto pode 
estudar-se por intermédio da permitividade complexa 
=" — je”, isto é, o deslocamento eléctrico não 
se encontra em fase com o campo eléctrico mas sim 
em atraso devido à inércia da polarização, que 
determina uma capacidade complexa C=(—j Ei. 
C, , sendo nestas condições [= jo(—j*)C, V 
e portanto com 1 =1, + ijl« deduz-se le = 
= 08" C Ve Je = mec, V donde o factor de 
perdas de polarização tg %p = Tait resulta tg 9,, = 


o, 
E 
! 


A presença da corrente activa corresponde de 
facto a uma condutividade dieléctrica de polarização 0 


f 
.— Ig? 
V-e 


p 
que se pode deduzir da lei d' Ohm Ip =G, V, 


sendo Gp a correspondente condutância Gp — 
=p A/d e juntamente com C, =: A/d se 
permite concluir 5; =wº, “ etg)p=op/ms É. 

Note-se porém que =” apenas tem valores apre- 
ciáveis, e portanto as perdas de polarização, 
quando a inércia da polarização for significativa, 
o que fisicamente equivale a um não acompanha- 
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mento do ritmo do campo aplicado, necessária- 
mente para elevadas frequências. Este fenómeno 
será objecto de análise detalhada. 


4. DIELÉCTRICO REAL 


Os materiais dieléctricos são porém mais ou 
menos imperfeitos, possuindo uma condutividade 
eléctrica não nula, embora muito reduzida, que 
provoca as denominadas perdas dieléctricas de con- 
dução iónica. 

Clássicamente o problema das perdas dieléc- 
tricas é encarado através do aparecimento de uma 
componente activa da corrente total pelo que, 
para v = Vsin q t, serái= Isin (»t + 2) onde 
o <=/2 radel=1p +j Iç com tg 9 = Iç/le- 
Assim a potência activa, relativa às perdas die- 
léctricas, deixa de ser nula, P = (1/T). É jp vidt= 
=V Icos 7/2 ou em valores eficazes, P = Vlcos 7 

Em geral caracterizam-se as perdas dieléctricas 
através do úngulo de perdas 3 == =/2 — v, donde 
P=VIsinó, com Icosy=Ip=-Ictgie 
lo =» CV, concluindo-se a expressão clássica 


P=oCVtgj (1) 


que traduz a proporcionalidade das perdas die- 
léctricas com a frequência w», a capacidade C do 
dieléctrico (ou a sua permitividade :), o quadrado 
da tensão aplicada e o factor de perdas tg 2. 

Este resultado é importante em alta tensão 
pois sendo V muito elevada, e para uma dada 
frequência (50 Hz no caso das redes de distri- 
buição de energia eléctrica) a redução das perdas 
dieléctricas só se pode conseguir pelas caracte- 
rísticas do material dieléctrico: por um lado pe- 
quena permitividade «: e por outro lado baixa 
condutividade (ou seja reduzido ângulo de per- 
das 9). Trata-se pois de um problema fundamen- 
tal da tecnologia dos materiais isolantes. 

Note-se a aproximação P = Ultgô no caso de 
grandezas sinusoidais, válida para pequenos valo- 
res de à pois então tg)=sin 3==cosy ou seja 
=. 

O esquema eléctrico equivalente do condensador 
pode conceber-se pelo paralelo de uma resistên- 
cia R de perdas de condução com a capacidade 
C, normalmente utilizado, mas nalguns casos 
interessa considerar o esquema série da resistên- 


cia R, com, havendo a relação RR, =1/(»C) ' 
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entre R e R, admitindo a igualdade das capaci- 
dades de ambos os esquemas como fácilmente se 
calcula pela equivalência das impedâncias dos 
dois esquemas. 

Atendendo, contudo, às perdas de polarização 
o esquema equivalente refere-se ao paralelo de 
R com a capacidade complexa C pelo que I = 
= VR tjw(:'—je”) C. V sendo a componente 


óhmica da corrente global Fi, = (1/R + me” C)V 
:&* C,V, donde 


=— (5) E 


e a componente reactiva E 
tgi= —=— + - (2) 


que se pode considerar tg) =tg ic +tg?p,a pri- 
meira parcela relativa às perdas de condução 
iónica tg)ç=1/»:'RC, e a segunda parcela 
correspondente às perdas de polarização tg 2 p = 
= 6“ lt, 


5. CONDUTIVIDADE DIELÉCTRICA 


Após o instante de aplicação de uma tensão 
continua a um dieléctrico constata-se em função 
do tempo uma variação da condutividade em 
quatro zonas (fig. 1-a), a saber: 


Zona I—Em menos de 107 * seg após o ins- 
tante de aplicação da tensão ao dieléctrico a 
condutividade decresce até se estabilizar devido 
à orientação dos dipolos da estrutura do mate- 
rial, definindo-se a condutividade de polarização. 


Zona Il —Entre 107” e 10 seg os electrões 
livres carregam os eléctrodos do condensador, 
correspondendo de facto à condutividade a 50 H.. 


Zona III — Para tempos de 10 a 10º seg dá-se 
uma forte redução da condutividade devida a 
duas causas: empobrecimento de rápidas partí- 
culas livres transportadoras de carga e grande 
concentração de cargas eléctricas nos eléctrodos. 
Poder-se-á referenciar este intervalo por conduti- 
vidade dieléctrica transitória. 


Zona IV —- Quando o tempo de aplicação da 
tensão contínua excede 10º seg atinge-se um 
valor constante que determina a chamada condu- 
tividade estacionária, 
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92 pr? 1 19º 


JA 


Cc) 


Fig. 1 — Condutividade dieléctrica, 
a) Em função do tempo. 
hj Em função da temperatura (escala linear), 
c) Em função do inverso da temperatura abso- 
luta (escala semi-logarítmica). 


A dependência da temperatura da condutivi- 
dade dieléctrica em corrente contínua e em 
corrente alternada à frequência industrial é re- 
gida pela lei de Van't Hoff 


. W 
“q = es “do exp (> = (3) 


onde 7,, é uma constante dependente da concen- 
tração iónica, W a energia térmica de activação, 
K a constante universal de Boltzmann e Ta tem- 
peratura absoluta, cuja representação gráfica 
(fig. 1-b, c) em escala semi-logarítmica em fun- 
ção de 1/T é uma recta pois ln 2, = In 0, — 
— W/KT. 


6. POLARIZABILIDADE 


De um modo genérico existem os seguintes 
mecanismos de polarização : 

1 — Polarização de orientação: Nos dieléctricos só- 
lidos e líquidos podem existir moléculas polares, 
isto é, cujos centros de gravidade das cargas eléc- 
tricas positivas e negativas não coincidem (exemplo: 
no ácido clorídrico existe numa deslocação de 
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carga eléctrica entre os átomos H e Cl). Sob o 
efeito de um campo eléctrico aplicado exterior- 
mente, estas moléculas rodam em torno dos seus 
eixo de simetria de modo a que cada dipolo 
tende a alinhar-se com o campo. 


2 — Polarização de distorsão: Nas moléculas pola- 
res o campo eléctrico, além do efeito de rotação 
que define a polarização de orientação, pode des- 
locar os átomos entre si alterando as distâncias 
que os separam, e por consequência o momento 
dipolar da molécula, fenómeno que constitui 
a polarização atómica. Outro tipo de polarização de 
distorsão pode verificar-se ainda nos átomos pois 
sendo constituídos por uma carga eléctrica posi- 
tiva central (núcleo) envolvida por uma núvem de 
electrões com determinadas simetrias definidas pe- 
los seus estados quânticos, quando um campo eléc- 
trico é aplicado os electrões são ligeiramente deslo- 
cados nas suas órbitas em relação aos núcleos 
induzindo um momento dipolar nos átomos 
que caracteriza a polarização electrónica, também 
chamada polarização óptica porquanto os dipolos 
produzidos pela deslocação elástica dos electrões 
têm frequências naturais iguais ou superior às 
da luz visível. Sob o ponto de vista técnico este 
tipo de polarização é de menor importância. 


3 — Polarização interfacial; Um dieléctrico sólido 
não é homogéneo em toda a sua extensão, 
havendo zonas cristalinas, zonas amorfas, impu- 
rezas, deslocações, etc. Partículas livres transpor- 
tadoras de carga eléctrica sob a acção de um campo 
eléctrico aplicado migram através da estrutura 
cristalina podendo acumular-se em defeitos estru- 
turais que com uma carga eléctrica imagem indu- 
zida num eléctrodo originam um momento dipolar. 
Nos dieléctricos líquidos com comportamento 
newtoneano não se verifica este tipo de polarização. 

O momento elementar de um dipolo consti- 
tuído por duas cargas eléctricas de diferente 
polaridade + Q à distância d é m== Qd e tem 
a direcção da carga positiva para a negativa. 
Pode no entanto assumir-se que este momento 
é proporcional a um campo eléctrico local interno 
ao dieléctrico E, através de uma constante de 
proporcionalidade x, chamada polarizabilidade do 


dipolo, sendo m = « E, e tendo 2 o significado 
do momento dipolar médio por unidade de 
intensidade de campo eléctrico. 

Considerando n dipolos elementares por unidade 
de volume com o momento elementar médio m 
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— > -— + 
o vector polarização será P == nm e daí a equação 


de Clausius P= n aE,. sendo a polarizabilidade « 
igual à soma das polarizabilidades individuais 
devidas a cada tipo de polarização: «= a, + 


à a + Ea Ei A. 
7. PERMITIVIDADE 


Um dipolo eléctrico que pode rodar livremente 
alinha-se com um campo aplicado e se este 
é alternado o dipolo roda ao ritmo do campo, 
ocorrendo uma ressonância quando a frequência 
do campo aplicado for igual à frequência de 
rotação natural do dipolo. Se porém o dipolo 
está num meio viscoso, provocando um certo 
amortecimento por atrito, a sua resposta ao 
campo eléctrico é do tipo relaxação, com uma 
certa constante de tempo, pelo que resulta um 
atraso de fase da rotação relativamente ao campo 
aplicado. 

Mas um dipolo eléctrico pode vibrar além 
disso pela variação da distância entre as cargas 
eléctricas devido ao campo aplicado, vibração 
encarada como um oscilador harmónico linear, 
dando-se também neste caso uma possível resso- 
nância ou relaxação de acordo com o amorteci- 
mento. 

Esta análise evidencia a dependência do valor 
da permitividade com a frequência, problema 
resolvido pela teoria de Debye. Além disso 
também depende da temperatura, esclarecendo-se 
igualmente pela mesma teoria o respectivo fenó- 
meno de relaxação. 


7.1 Influência da frequência na permitividade 


a) Devido à polarização de orientação 

Quando um dieléctrico é submetido a um 
campo eléctrico alternado o fenómeno da polari- 
zação apresenta uma certa inércia que origina 
um certo atraso do deslocamento eléctrico em 
relação ao campo eléctrico aplicado, o que se 


traduz por uma permitividade complexa: D = 


= e (é —j:”) E. Para elevadas frequências a 
orientação dos dipolos ao ritmo do campo não se 
pode realizar devido à inércia, pelo que :” se anula. 

Este comportamento em função da frequência 
exprime-se analiticamente por uma relação do 
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tipo é (») =: + [Pa (t). exp (jwt) dt, onde 


6) 

E. é o valor da permitividade à frequência infinita, 
sendo ôbviamente constante, e a (t) o factor de 
queda de : (») em consequência do atraso da 
polarização em referência ao campo, o que 
equivale à redução da polarização quando se retira 
repentinamente um campo estacionário aplicado 
ao dieléctrico ou à subida da polarização quando 
se aplica repentinamente esse campo estacionário, 
portanto num regime transitório caracterizado 
por um tempo Z chamado tempo de relaxação. 

Debye propôs uma forma exponencial para 
este factor de queda da permitividade «a (t) = 
== « (0). exp (— t/Z), onde o tempo de relaxação 
Z, característico do dieléctrico, pode ser função 
da temperatura mas não do tempo. No caso da 
corrente contínua, w==0, seja e (0) =: e por 
outro lado : =: +z.«(0) donde a (t) = 
= [(6,—€,.)/7]. exp (-t/0). 

Integrando c(m) == :.., + Ja [(e, — too)/7]. 
«exp [(j wo — 1/Z) t] dt obtém-se : (w) = “so + 
+ (1, — Eoo)/(1 — jmT)e, por separação das com- 
ponentes real e imaginária, as chamadas equações 
de Debye 


E — É 
== to + +, (4) 
1 + (mz) 


(e, ca oo) 6) (é 


1 + (mz) 


Ef mam 


(5) 


Pela representação gráfica (fig. 2-a) de “& (mm) 
e <” (w) verifica-se a constância de : para baixas 
frequências (Zona 1) e para muito altas frequên- 
cias (Zona III) havendo uma zona intermédia 
de relaxação na qual se manifesta «” (Zona Il). 


b) Devido à polarização de distorsão 


A polarização de distorsão resulta da desloca- 
ção de electrões (polarização óptica) e de átomos 
(polarização atómica) por acção de um campo 
eléctrico aplicado. 

Basicamente as cargas eléctricas elasticamente 
ligadas possuem uma posição de equilíbrio que será 
retomada quando o campo aplicado deixar de exis- 
tir por meio de oscilações amortecidas em torno 
da correspondente posição de equilíbrio. Este 
modelo permite antever uma ressonância, na 


313 


irc 
H 
E (cu) | | 
LA — ,.. meme ... — — e gre 
2 Bl. intarfe fot. orient, ol. dust, tu 
4) 


Fig. 2 — Variação da permitividade de um dieléctrico 
com a frequência, 
a) Devido à polarização de orientação. 
b) Devido à polarização de distorsão., 
c) Esquema de um dieléctrico heterogéneo, 
dy Variações de :' e «” com a frequência, 


qual se absorve uma energia máxima, se o campo 
corresponde à frequência própria. 

Considerando portanto o comportamento das 
cargas eléctricas como osciladores harmónicos 
lineares a equação do movimento na direcção x 
será x“ +2A mw x'+uwi x=0, com o factor 
de amortecimento À e a frequência angular natu- 


ral », , cuja solução, para pequenos amorteci- 
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mentos, é do tipo x=x,. exp (— At) . cosm,t, 
definindo-se o tempo de relaxação Z do oscila- 
dor por T=-1/A correspondente ao tempo em 
que a amplitude de oscilação se reduz a 1/e do 
seu valor máximo. 

Frcehlich propôs uma redução da polarização 
semelhante à de Debye afectada de cos (mt + “), 
onde ! é o ângulo de fase relativo ao atraso da 
polarização em relação ao campo aplicado, que 
permite concluir sobre a ressonância (fig. 2-b): 
apresenta-se um máximo *, =“. + (A:/4).m,7 
e um mínimo com (dE/4) Tt, com a 
diferença As entre os valores de :' antese depois 
do período de relaxação, enquanto “qax = 
= (A:/2).mw z para a frequência », de modo 


aque: /2 corresponde às frequências», T1/Z. 


! 


c) Devido à polarização interfacial 


A heterogeneidade de um dieléctrico refere-se 
na sua forma mais simples (num caso mais 
complexo o raciocínio não se altera em princípio) 
a duas camadas, com uma superfície de discon- 
tinuidade, caracterizadas pelas permitividades 
:, e &,  condutividades 7,e7,eespessurasd ed,, 
respectivamente. 

Quando um campo é aplicado o comporta- 
mento do dieléctrico traduz-se pela admitância 
equivalente Y=Y,..Y/Y, +HY)=jr(” — 
—j:*)C;,,comt,=R,C,,7%=-R,C, et = 
=(R,2,+R,Z,)/R,sendoR=R, +R,, que 
conduz às relações seguintes: 


o — E 


E E É as + -— , —s (6) 
1 + (mz) 


que é exactamente igual à expressão de : no caso 
da polarização de orientação (4), quer dizer, 
a medida da parte real da permitividade relativa 
dá a soma das contribuições das polarizações 
interfacial e de orientação, não sendo possível 
separá-las ; 

1 (: 


E AR ad 2. E o ) mt (7) 
mRC, 1 + (nz) 


que apresenta um termo adicional ao valor de 
=” no caso da relaxação de Debye (5), inversa- 
mente proporcional a » (tende para infinito quando 


m tende para zero numa evolução hiperbólica), 
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isto é, a polarização interfacial pode distin- 
guir-se da relaxação por orientação pela obser- 
vação da variação de «” abaixo da frequência de 
relaxação. 


7.2 Influência da temperatura na permiti- 
vidade 


Num dieléctrico líquido (caso prático: óleo de 
transformadores) não sujeito a nenhum campo 
eléctrico os dipolos têm uma distribuição tal que 
o momento resultante é nulo. Se se aplica numa 
dada direcção um campo eléctrico constante então 
os dipolos elementares orientam-se de maneira a 
que o momento total deixa de ser nulo. Se se 
voltar a anular o campo eléctrico, em virtude do 
movimento térmico molecular (Browniano) o 
momento global torna outra vez a anular-se. 

Sob as hipóteses de Debye das moléculas serem 
esféricas de raio a, de um só tipo, sem acção entre 
sie num fluido viscoso contínuo não polar de visco- 
sidade 7», através da lei de Stockes, que define a 
constante de atrito B do movimento das molé- 
culas dipolares (sendo o momento m produzido 
por um campo eléctrico fazendo um ângulo ? com 
o eixo do dipolo de expressão m = B d 9/dt) com o 
valor B==87 na, deduz-se para o tempo de 
relaxação, tendo em atenção ao movimento 
Browniano que afecta inevitavelmente o ritmo de 
variação do ângulo com o campo, TZ = B/2k T, 
ou seja 


5 
axna 

t= (8) 
k T 

expressão que indica o resultado interessante de 

que a variação do tempo de relaxação Z com a 

temperatura é proporcional à relação do coeficiente 

de viscosidade e da temperatura absoluta. 
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Fig. 3 — Variação da permitividade de um dieléctrico 
com a temperatura. 
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À representação gráfica (fig. 3) de é e :” em 
função da temperatura apresenta três zonas 
distintas : 


Zona I— para baixas temperaturas a viscosi- 
dade é muito elevada pelo que os dipolos estão 
«congelados» ; 

Zona II— para temperaturas intermédias a 
viscosidade permite a orientação dos dipolos mas 


—+ —+> 
o campo eléctrico E não está em fase com D pelo 


que surgem perdas de polarização (W = E D/2) : 
Zona III — para altas temperaturas o coeficiente 


de viscosidade é muito reduzido pelo que E está 


em fase com D, sendo tg 9p =. 


8. PERDAS DIELÉCTRICAS 


Em corrente alternada podem ocorrer três tipos 
de perdas dieléctricas: perdas de condução iónica, 
relativas a uma autêntica condutividade eléctrica, 
perdas de polarização, devidas aos diferentes 
aspectos da polarização do dieléctrico sob a acção 
do campo eléctrico aplicado, e perdas de ioni- 
zação, para tensões suficientemente elevadas 
e causadas por descargas parciais. 


8.1 Perdas de condução 


No caso de líquidos e gases estas perdas resultam 
do movimento migratório permanente de electrões 
e iões. No caso de sólidos as partículas embora 
não emigrem sofrem deslocações. 

Para a mobilidade b das partículas transpor- 
tadoras da carga q e de concentração N a 
condutividade será 3, = N bq sendo a densidade 


de corrente J=c,E ou ]J==NbqE e as perdas 


específicas (por unidade de volume) pç = E ou 


Po=4 EP e as perdas totais de condução do 
dieléctrico de volume U de valor P; = /y pç dU. 

Para domínios de temperatura não muito 
grandes é válida a lei de Van't Hoff da condutivi- 
dade e portanto para determinados campos eléc- 
tricos (ou para uma tensão constante) as perdas 
variam com a lei exporencial de q,. 

De facto, como se deduziu, as perdas relativas 
à condução caracterizam-se pelo factor de perdas 
tg)c=1/»RC, ou tgiç=04,/mn:,€, quer 
dizer, as perdas de condução iónica são directa- 
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